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Abstract. Pinwheel patterns and their higher dimensional generalisations display continuous 
circular or spherical symmetries in spite of being perfectly ordered. The same symmetries 
VP ' show up in the corresponding diffraction images. Interestingly, they also arise from amorphous 

systems, and also from regular crystals when investigated by powder diffraction. We present 
C ~ 3 ■ first steps and results towards a general frame to investigate such systems, with emphasis on 

£*\} ' statistical properties that are helpful to understand and compare the diffraction images. We 

concentrate on properties that are accessible via an alternative substitution rule for the pinwheel 
(^J , tiling, based on two different prototiles. Due to striking similarities, we compare our results with 

1 the toy model for the powder diffraction of the square lattice. 

1. Pinwheel patterns 

CM 

The Conway-Radin pinwheel tiling 14 , a variant of which is shown in Figure[TJ is a substitution 
tiling with tiles occurring in infinitely many orientations. Consequently, it is not of finite local 
complexity (FLC) with respect to translations alone, though it is FLC with respect to Euclidean 
motions. This property distinguishes the pinwheel tiling from the majority of substitution tilings 
considered in the literature. As a consequence, its diffraction differs considerably from that of other 
tilings, and despite a growing interest in such structures [131 1121 HI I18 | . the diffraction properties 
have only been partially understood to date. 

Whereas the pinwheel tiling is the most commonly investigated example, there are other tilings 
with infinitely many orientations, compare [15} for an entire family of generalisations. Yet another 
example is shown in Figure [3 It has a single prototile, an equilateral triangle with side lengths 
1, 2 and 2. Under substitution, the prototile is mapped to nine copies, some rotated by an angle 
9 = arccos(l/4), which is incommensurate to ir (i.e., 8 ^ 7rQ). Thus, the corresponding rotation 
Re is of infinite order, and the tiles occur in infinitely many orientations in the infinite tiling. 
Here and below, R a denotes the rotation through the angle a about the origin. More examples of 
$H I tilings with tiles in infinitely many orientations can be found in [Zj. 

It was shown constructively in 12 that the autocorrelation 7 of the pinwheel tiling has full cir- 
cular symmetry, a result that was implicit in previous work (14] . As a consequence, the diffraction 
measure 7 of the pinwheel tiling shows full circular symmetry as well. To make this concrete, we 
now construct a Delone set from the tiling. Recall that a Delone set A in Euclidean space is a 
point set which is uniformly discrete (i.e., there is r > such that each ball of radius r contains 
at most one point of A) and relatively dense (i.e., there is R > such that each ball of radius R 
contains at least one point of A). Let T be the unique fixed point of the pinwheel substitution of 
Figure n that contains the triangle with vertices (§,— \,— i )> ( — §' !)■ This fixed point T is 
the same as the one considered in ^2] ■ We now define the set of control points At of T to be the 
set of all points u + + U ~ 4 W such that the triangle with vertices u, v, w is in T and uv is the 
edge of length one. This choice of control points is indicated in Figure ^ (left) and is the same as 
in P2]. 

Recall that the natural autocorrelation measure of a Delone set A is defined as 

x,y£AnB R 

where the limit is taken in the vague topology and exists in all examples discussed below; for 
details, see [21 El US]- Here, S x denotes the Dirac measure in x, and Br the closed ball of radius R 
centred at the origin. The Fourier transform 7 is then the diffraction measure of A, whose nature 
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Figure 1 . The pinwhccl substitution rule and a patch of the pinwheel tiling T. 
The points in the left part indicate how the point set A = Aq- arises from the 
pinwhccl tiling. 

is often the first property to be analysed. Since 7 is a translation bounded measure on R 2 , it has 
a unique decomposition, relative to Lebesgue measure, into three parts, 

1 Tpp 1 7sc ' lac? 

where the pure point part 7 pp is a countable sum of (weighted) Dirac measures, 7 ac is absolutely 
continuous with respect to Lebesgue measure, and 7 5C is supported on a set of Lebesgue measure 
0, but vanishes on single points. 

It was shown in J2| that the autocorrelation 7^ of the pinwheel control points A = Aq- satisfies 

(2) 1a = s o + Vi^fh- = v{r)fJ-r, 

rev\{o} rev 

where P is a discrete subset of [0, +00), fi r denotes the normalised uniform distribution on the 
circle rS 1 — {x £ R 2 \x\ = r}, and r/(r) is a positive number. Note that /i = S . In particular, j A 
shows perfect circular symmetry, as docs the diffraction measure "f A . This settles the pure point 
part: Since j A is a translation bounded measure, and the Fourier transform of such a measure 
is also translation bounded, it follows from the circular symmetry that there are no Bragg peaks 
except at 0. Moreover, a standard argument 9; gives 

7 PP = (dens(yl)) 2 (5 = 5 , 

because the density dens(/l), i.e., the average number of points of A per unit area, is 1. This 
follows from the fact that, in our setting, each triangle has unit area and carries precisely one 
control point. 

Below, we give more detailed information about T> and r](r), which is needed to shed some light 
on the nature of 7 SC and 7 ac . 

Proposition 1. The pinwheel Delone set A as defined in Figure^ satisfies: 

(i) A C \ \ neJj R n el? , where :— 2arctan(|). 

(ii) Ac{(f,f)|m,neZ, fc£N }. 

(iii) The distance set T) — V A := {\x — y\ \ x, y £ A} is a subset of {y p ^ t q \ p,q,£ G N } . 
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Figure 2 . Another substitution rule which generates a tiling with circularly sym- 
metric autocorrelation. The decoration of the triangles indicates that, in contrast 
to the classic pinwheel, no enantiomorphic pairs of triangles occur here. 



In plain words, A is a uniformly discrete subset of a countable union of rotated square lattices, 
all elements of A have rational coordinates, and, as a consequence, all squared distances between 
points in A are rational numbers of the form (p 2 + q 2 )/5 e . The fact that A is supported on such 
a simple set is an interesting property of the pinwheel tiling. It is not clear whether a similar 
property, for a suitable choice of control points, can be expected for other examples, such as that 
of Figure 

These results were obtained by means of an alternative substitution, the kite domino substi- 
tution shown in Figure |21 which generates the same Delone set A. The kite domino substitution 
is equivalent to the pinwheel substitution in the sense that the corresponding tilings are mutually 
locally derivable (MLD) in the sense of 0|, i.e., they can be obtained from each other by local 
replacement rules. Moreover, the Delone set A is MLD with both tilings. 

Because of the strong linkage between the diffraction spectrum of a Delone set A and the 
dynamical spectrum of the associated dynamical system (X(/l),K rf ), we consider the hull of A 

I DT 

X(A) = M? + A , 

where completion is with respect to the local rubber topology (LRT), see 3 and references therein 
for details. Roughly speaking, and restricted to the special case under consideration, this means 
that X(A) contains all translates of A and all Delone sets which are locally congruent to some 
translate of A. 

In addition to Proposition [T] the kite domino substitution gives access also to the frequency of 
configurations in the tilings. The frequency of a finite set L C A is defined as 

freq(X) = lim ——7 card{F C A n Br I F is congruent to L\. 

Note that this definition is up to congruence of the finite sets, not up to translation (which is not 
reasonable here). The frequency module of A is the Z-span of {freq(L) | L C A finite}. 

Proposition 2. The frequency module of A is { 26 " 5 <; | fn G Z, I S N }. It is the same for all 
elements ofK(A). In particular, all rj(r) in @ are rational. 

Using the kite domino substitution, one can determine some frequencies of small distances 
exactly. These are given below, together with some other values (marked by an asterisk) where 
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Figure 3. The kite domino substitution rule and a patch of a kite domino tiling. 
This patch is equivalent to the pinwheel patch in Figure ^ The dots in the patch 
indicate how A arises from the tiling. 

the frequencies are estimated by analysing large approximants of the pinwheel tiling. 
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By using 'collared' tiles (which refers to the 'border- forcing' property of |11|L one can in principle 
derive all frequencies in D, see In fact, the frequency of pairs of points with distance r = 1 
in the table above was calculated this way. However, the computation of each single frequency 
requires a considerable amount of work, and a closed formula for all frequencies seems out of reach. 

2. Diffraction 

The diffraction of a crystal which is supported on a point lattice in R d is obtained by the Poisson 
summation formula for Dirac combs [SHU- If J 1 is a lattice, the autocorrelation of the lattice Dirac 
comb Sr is dens(-T) <5r, and the diffraction measure reads 

(4) S r = (densOT)) 2 .* r ., 

where r* is the dual lattice of r. 

A radial analogue of Eq. Q is derived in . It is an analogue of the Hardy-Landau- Voronoi 
formula |10| in terms of tempered distributions. Let us first explain this for the example of the 
square lattice Z 2 . Let V D be the distance set of Z 2 , and r] a {r) :— card{x 6 Z 2 \ \x\ — r} the 
shelling numbers of Z 2 , see 2 for details. Then 

(5) Vn( r )Vr) = ^a( r )K= Vair)^ 

rev n reV u reV n 

with [i r as above. Again, the sum is to be understood as a vague limit. The fact that Z 2 is 
self-dual as a lattice (i.e., (Z 2 )* = Z 2 ) implies that the same distance set enters all three sums in 
(|5j). In the general case, with an arbitrary lattice r, one has to use the distance set of the dual 
lattice r* of r, in analogy with Q). This gives the following result PP, valid in Euclidean space 
of arbitrary dimension. 
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Theorem 1. Let r be a lattice of full rank in R d , with dual lattice r* . If the sets of radii for 
non-empty shells are T> r and T> r ,, with shelling numbers rj r (r) — cardjx G -T | \x\ = r} and 
r] r , (r) defined analogously, the classical Poisson summation formula has the radial analogue 

(6) Vr( r )f-r) = dens(r) Y Vr*(r)f- r , 

rev r rev r , 

where [i r denotes the uniform probability measure on the sphere of radius r around the origin. 

Eq. © also implies 

( 7 ) (Y IrW/fr) = Vr( r )K 

rev r rev r 

in the sense of tempered distributions. This equation is useful in numerical calculations of pinwhcel 
diffraction spectra. 

2.1. Pinwheel diffraction. Let us return to the pinwheel pattern A. In view of Eq. (J2J) in 

connection with Proposition ^ the sum in Eq. (J2J) can be recast into a double sum, 

oo 

(8) 7 = Y Y %( r )Mr> 

where the choice of the T)g(r) is not unique, but restricted by the condition X)fco r le( r ) = r l{ r )- If 
the T) t (r), for fixed £, could be chosen to be 'lattice-like' — in the sense that they form a sequence 
of shelling numbers for some lattice — we were in the position to apply Eq. © to each inner sum 
in Eq. © individually. Observing that 

(9) (5-^ 2 Z 2 )* = h 1 ' 2 Z 2 

this would give rise to terms of the form ^2 re5 e/2 V Ve( r ) V-r- D ue to the continuity of the Fourier 
transform on the space of tempered distributions, this would imply the diffraction of the pinwheel 
tiling to be purely singular. 

However, things are not that simple in the case of the pinwheel tiling. In particular, the rj^fr) 
cannot be chosen to be lattice-like, as a consequence of the Delone nature of A. Nevertheless, it 
may still be possible to find coefficients ei, not necessarily positive, such that the autocorrelation 
of A can be written as 

T = Y e e Y ^W^r- 

This general form permits continuous parts in the diffraction different from the ones arising from 
©: There may be singular continuous parts apart from {rS 1 | r G T> \ {0}}, and even absolutely 
continuous parts. In fact, numerical computations indicate £Q the presence of an absolutely 
continuous part in the diffraction. 

Independent of an affirmative answer of the open questions, there is a striking resemblance with 
the powder diffraction image of the square lattice. This is a consequence of Proposition^] So let 
us close this article with a simplified approach to powder diffraction patterns, and a comparison 
of the two images. 

2.2. Square lattice powder diffraction. Instead of performing a diffraction experiment with 
a large single crystal, it is often easier to use a probe that comprises many grains in — ideally — 
random and mutually uncorrelated orientations |17| . This setting can be modelled mathematically 
as follows. Let R be a rotation of infinite order, i.e., R n ^ id for all integer n ^ 0. It follows from 
Weyl's lemma that {R n x n G N } is uniformly distributed on the unit circle S 1 for any x G S 1 . 
For simplicity, we require Rl? (~1 Z 2 — {0}. Then, for large N, the set 

i N 
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Figure 4. Numerical approximation to the radial intensity structure I(k) of 
the pinwheel diffraction (solid line) in comparison with the powder diffraction 
structure of the square lattice (grey bars). The relative scale has been chosen 
such that the heights of the first peaks at k = 1 match. Note that the central 
intensity is suppressed. 



can serve as an idealised powder emerging from a two-dimensional crystal supported on Z 2 . Note 
that the prefactor 1/N appears because we idealise the arrangement of disoriented grains as an 
overlay of mutually rotated infinite copies of Z 2 . 

It is not hard to show [J that the corresponding autocorrelation is given by 




where A denotes Lebesgue measure in R 2 . The limit of the bracketed term, as N — > oo, shows 
perfect circular symmetry, and is a reasonable approximation of the powder autocorrelation. An 
application of Theorem \\\or Eq. (0), yields 



N-l 



(10) n^(n <W 2 ) = (^2 Vn( r ) Vr) = ^2vn( r )Vr 



j=0 rev. 



a 



This shows that, for an ideal square lattice powder as above, one can expect a diffraction image 
that, beyond the central intensity, consists of concentric rings of radius r 6 T> D with total intensity 
Tj n (r). In this simplified version, there is no absolutely continuous part, though this would be 
present in a more realistic model. 

To compare the powder diffraction of Z 2 with the pinwheel diffraction, we note that it is 
sufficient to display the radial structure. Furthermore, we cannot compare the central intensity, 
wherefore we suppress it in both cases. While Eq. I|1(J|I gives a closed formula for the intensity 
of the rings in the powder diffraction, we currently only have a numerical approximation to the 
pinwheel diffraction, based on Eq. |J7J and a large patch, see £Q for details. Figure which 
speaks for itself, shows striking similarities of the singular parts. It is thus plausible that further 
investigations in this direction might ultimately reveal the full nature of the pinwheel diffraction. 
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